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METHOD AND APPARATUS FOR ASSAYING A DRUG CANDIDATE TO 
ESTIMATE A PHARMACOKINETIC PARAMETER ASSOCIATED THEREWITH 



5 assaying a drug candidate and, more specifically, to a method for measuring the binding 
interaction between a drug candidate and sensing surface-bound biomolecules of a 
biosensor to determine a binding interaction parameter of the drug candidate, and then 
comparing the binding interaction parameter against a predetermined drug correlation 
graph {e.g., a mathematical expression) to estimate at least one pharmacokinetic 
1 0 parameter. 



chemical, physical and biological properties associated with low molecular weight 
substances, particularly in the context of drug discovery. For example, researchers are 
15 often concerned with determining a variety of chemical, physical and biological 
properties associated with drug candidates for screening purposes. The determination 
of such properties often plays a pivotal role in the drug development and screening 
process. 



20 duration of the pharmacological effect of a systemically acting drug are functions not 
only of the intrinsic activity of the drug, but also of its absorption, distribution, 
metabolism, and excretion (ADME) characteristics within the human body. These so- 
called ADME characteristics are all intimately related to the scientific discipline known 
as "pharmacokinetics." Pharmacokinetics is commonly referred to as the study of the 

25 time courses (i.e., kinetics) associated with the dynamic processes of ADME of a drug 
and/or its metabolites within a living organism, and is closely interrelated with the 
fields of biopharmaceutics, pharmacology, and therapeutics. 




This invention is generally directed to a method and apparatus for 



BACKGROUND OF THE INVENTION 



A variety of experimental techniques are currently used to determine 



More specifically, it has long been recognized that the intensity and 



Because the body delays the transport of drug molecules across 
membranes, dilutes them into various compartments of distribution, transforms them 
into metabolites, and eventually excretes them, it is often difficult to accurately predict 
the pharmacological effect of promising new drug candidates. Researchers, however, 
5 commonly use pharmacokinetic ADME studies as one method to predict the efficacy of 
a drug at a site of action within the body. 

Traditionally, researchers involved with preclinical ADME studies have 
used pharmacokinetic/mathematical models coupled with actual drug concentration data 
from blood (or serum or plasma) and/or urine, as well as concentration data from 

10 various tissues, to characterize the behavior and "fate" of a drug within living 
organisms. As is appreciated by those skilled in the art, the mathematical equations 
associated with pharmacokinetics are generally based on models that conceive the body 
as a multicompartmental organism. In such models it is presumed that the drug and/or 
its metabolites are equitably dispersed in one or several fluids/tissues of the organism. 

15 Any conglomerate of fluid or tissue which acts as if it is kinetically homogeneous may 
be termed a "compartment." Each compartment acts as an isotropic fluid in which the 
molecules of drug that enter are homogeneously dispersed and where kinetic 
dependencies of the dynamic pharmacokinetic processes may be formulated as 
functions of the amounts or concentrations of drug and metabolites therein. Stated 

20 somewhat differently, the conceptual compartments of the body are separated by 
barriers that prevent the free diffusion of drug among them; the barriers are kinetically 
definable in that the rate of transport of drug or metabolite across membrane barriers 
between compartments is a function of, for example, the amounts or concentrations of 
drug and metabolites in the compartments, the permeability of various membranes, 

25 and/or the amount of plasma protein binding and general tissue binding. 

More specifically, pharmacokinetic/mathematical models are commonly 
used by pharmacokineticists to represent drug absorption, distribution, metabolism, and 
excretion as functions of time within the various tissues and organs of the body. In such 
models, the movement of the administered drug throughout the body is concisely 

30 described in mathematical terms (e.g., a set of differential equations). The predictive 
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capability of such models lies in the proper selection and development of mathematical 
functions that parameterize the essential factors governing the kinetic process under 
consideration. 

For example, a drug that is administered by intravenous injection may be 
5 assumed to distribute rapidly in the bloodstream. A pharmacokinetic/mathematical 
model that describes this situation may be a tank containing a volume of fluid that is 
rapidly equilibrated with the drug. Because a fraction of the drug in the body is 
continually eliminated as a function of time (e.g., excreted by the kidneys and 
metabolized by the liver), the concentration of the drug in the hypothetical tank may be 
10 characterized by two parameters: (1) the volume of fluid in the tank that will dilute the 
3 drug, and (2) the elimination rate of the drug per unit of time, both of which are 

generally considered to be constant. Thus, if a known set of drug concentrations in the 
tank is determined at various time intervals by, for example, sampling, then the volume 
of fluid in the tank and rate of drug elimination may be estimated. This information 
15 may then, in turn, be used for predicting the disposition of the drug within a human 
body. 

Theoretically, an unlimited number of models may be constructed to 
describe the kinetic processes of drug absorption, distribution, metabolism, and 
excretion within the various tissues and organs of the human body. In general, 
20 however, the number of useful models is limited due to practical considerations 
associated with blood, tissue and/or organ sampling. As a result, and as is appreciated 
by those skilled in the art, two major types of models have been developed by 
pharmacokineticists: (1) compartmental models; and (2) physiologic models. 

In pharmacokinetic compartmental models, the body is represented as a 
25 series of compartments that communicate reversibly with each other. Each 
compartment is not a real physiological or anatomic region; rather, each compartment is 
considered to be inclusive of all tissues that have similar blood flow and drug affinity. 
For example, a compartmental model may consist of one or more peripheral 
compartments representing tissue(s) connected to a central compartment representing 
30 the blood stream. Conceptually, the drug moves dynamically into and out of the central 



4 



2 - s 



compartment and into and out of each of the peripheral compartments. As such, rate 
constants may be used to represent the overall rate process for the drug's disposition 
within each compartment. Compartment models are generally based on linear 
assumptions using linear differential equations, and are particularly useful when there is 
5 little information known about the tissues and their respective drug concentrations. 

In contrast, pharmacokinetic physiologic models are based on known 
anatomic and physiologic data, data which is kinetically described in view of the actual 
blood flow volumes responsible for distributing the drug to the various parts of the 
body. Because there are many tissue organs in the body, each tissue volume must be 

10 estimated and its drug concentration and rate of change described mathematically 
(tissues having similar blood perfusion properties, however, are typically grouped 
together). Unfortunately, much of the information required to adequately describe such 
pharmacokinetic physiologic models are often very difficult to obtain experimentally. 
Nevertheless, such physiologically based models are commonly used in conjunction 

15 with animal data and interspecies scaling techniques to predict the drug's disposition 
within a human body. 

More importantly, however, is that pharmacokinetic/mathematical 
models, and knowledge of their associated ADME parameters play an extremely 
important role in drug discovery and development. A typical example is a drug that is 

20 active following intravenous administration but is considerably less active after 
comparable oral doses. Having appropriate pharmacokinetic information may reveal (1) 
whether the drug was poorly absorbed to yield subtherapeutic circulating levels, or (2) 
whether the drug experienced presystemic metabolism to an inactive metabolite. Such 
information may also provide guidance for subsequent decisions, such as (1) whether to 

25 improve drug absorption by altering the salt form or formulation, (2) whether to 
investigate the possibility of making prodrugs, or (3) whether to consider a parenteral 
route of administration. 

In addition to the foregoing, pharmacokinetic/mathematical models are 
also generally considered extremely useful for, among other things: (1) predicting 

30 plasma, tissue, and urine drug levels with any dosage regimen; (2) calculating the 
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optimum dosage regimen for an individual patient; (3) estimating the possible 
accumulation of drugs and/or metabolites; (4) correlating drug concentrations with 
pharmacologic and toxicologic activity (i.e., pharmacodynamics); (5) evaluating 
differences in the rate or extent of availability between formulations (i.e., 
5 bioequivalence); (6) describing how changes in physiology or disease affect the 
absorption, distribution, and/or elimination of the drug; and (7) explaining drug-drug 
and food-drug interactions. 

Lastly, pharmacokinetic ADME data has also become an integral part of 
the pharmacological characterization process of promising new drug candidates. 
10 Regulatory agencies, such as the U.S. Food and Drug Administration, now require (1) a 
determination of pharmacokinetic ADME data in Phase I drug studies, and (2) a 
submission of pharmacokinetic ADME data as part of a New Drug Application. In this 
context, such pharmacokinetic ADME data is deemed essential for predicting the 
behavior and fate of the drug candidate within the human body, 
a ' 15 Accordingly, there is a need in the art for improved methods for 

Iq determining one or more pharmacokinetic parameters associated with absorption, 

^ distribution, metabolism, and excretion of a drug candidate. There is also a need for 

Q apparatuses useful for carrying out such methods. The present invention fulfills these 

needs and provides further related advantages. 

20 SUMMARY OF THE INVENTION 

In brief, the present invention is directed to a method and apparatus for 
assaying a drug candidate. More specifically, this invention discloses a method for 
measuring the binding interaction between a drug candidate and sensing surface-bound 
biomolecules of a biosensor to determine a binding interaction parameter of the drug 

25 candidate, and then comparing the binding interaction parameter against a 
predetermined drug correlation graph to estimate at least one pharmacokinetic 
parameter of the drug candidate. The at least one pharmacokinetic parameter may, for 
example, be one or more of ADME. 



In another embodiment of the present invention, at least two 
pharmacokinetic parameters of the drug candidate are determined, and in yet another 
embodiment, at least one pharmacokinetic parameter and a solubility property of the 
drug candidate are determined. Such pharmacokinetic parameters and/or solubility 
5 property may be determined when the one or more sensing surface-bound biomolecules 
are selected from, for example, liposomes, plasma proteins, CYP 450 enzymes, 
metabolic enzymes, or transport proteins. 

The biosensor used in the practice of the present invention may utilize a 
mass-sensing technique, such as surface plasmon resonance. In addition, the biosensor 
f «l 10 may further employ a sensor chip, wherein the sensor chip comprises a free electron 
\5 metal that includes a sensor surface, wherein the free electron metal is copper, silver, 

fU aluminum or gold. The sensor chip may further comprise a hydrogel coupled to the 

sensor surface, wherein the hydrogel has a plurality of functional groups, and wherein 
the one or more sensing surface-bound biomolecules are covalently bonded to the 
15 hydrogel. 

In a more specific embodiment, a sensor surface adopted for use with a 
biosensor is disclosed. The sensor surface comprises a hydrogel matrix coating coupled 
to a top surface of the sensor surface, wherein the hydrogel matrix coating has plurality 
of functional groups. At least two different liposomes are bonded to the plurality of 
20 functional groups at discrete and noncontiguous locations on the hydrogel mixtrix 
coating of the sensor surface. In one embodiment, the sensor surface is a sensor chip, 
and a free electron metal is interposed between the hydrogel matrix and the top surface 
of the sensor surface. 

In another embodiment of this invention, an apparatus is disclosed for 
25 assaying a drug candidate, wherein the apparatus comprises a biosensor having one or 
more sensing surface-bound biomolecules associated therewith and capable of 
measuring at least one binding interaction parameter of the drug candidate, and a 
computer memory containing a data structure for comparing the at least one binding 
interaction parameter against at least one mathematical expression correlated from 



m 



ill. 



• # 

7 

binding interaction data associated with known drug compounds to determine an 
estimate of at least one pharmacokinetic parameter of the drug candidate. 

In yet a further embodiment, a computer memory containing a data 
structure useful for assaying a drug candidate in accordance with the methods of the 
present invention is disclosed (as well as a generated data signal conveying the same). 
The data structure may be used to determine an estimate of at least pharmacokinetic 
parameter of the drug candidate. 

These and other aspects of the present invention will be evident upon 
reference to the following detailed description and related Figures. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1A illustrates a sensorgram reflective of the steady-state binding 
levels associated with a selected drug-biomolecule interaction. 

Figure IB illustrates an enlarged portion of the sensorgram of Figure 1A 
and shows sensorgram irregularities. 

Figure 1C illustrates an enlarged portion of the sensorgram of Figure 1A, 
wherein bulk-refractive index effects have been eliminated. 

Figure 2A illustrates a single "dip" depicting the reflected light intensity 
associated with a homogeneous sensor surface. 

Figure 2B illustrates a number of "dips" depicting the non-averaged 
reflected light intensities associated with a non-homogeneous sensor surface. 

Figure 2C illustrates a broadening of the "dip" depicting the averaged 
reflected light intensities associated with a non-homogeneous sensor surface. 

Figure 3 shows a high-level block diagram of an exemplary computer 
system for assaying a drug candidate in accordance with the methods of the present 
invention. 

Figure 4 depicts a correlation graph having known KD values for known 
compounds plotted along the abscissa (i.e., the x-axis) and corresponding measured KD 
values obtained via the biosensor plotted along the ordinate (i.e., the y-axis) for nine 
drugs with known levels of plasma protein binding 
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Figure 5 depicts a correlation graph having respective binding levels at 
100 |iM (R 100 jaM) divided by the molecular weight of known drug compounds' 
plotted along the abscissa (i.e., the x-axis) and corresponding human serum albumin 
binding percentage, as measured by equilibrium dialysis, plotted along the ordinate (i.e., 
5 the y-axis). 

Figure 6 depicts reference subtracted sensorgram traces for each of three 
drug candidates A-C. 

Figure 7 depicts a correlation graph having known fraction absorbed in 
humans (FA%) plotted along the ordinate (i.e., the y-axis) and corresponding calibrated 
10 (i.e., reference subtracted) steady state binding levels for each drug at 500 |iM plotted 



*3 along the abscissa in a 10-logarithm scale (i.e., the x-axis). 

DETAILED DESCRIPTION OF THE INVENTION 
>0 The present invention is directed to a method for assaying a drug 

a candidate and, more specifically, to a method and apparatus for measuring the binding 

1 5 interaction between at least one drug candidate and sensing surface-bound biomolecules 
O of a biosensor to determine a binding interaction parameter of the drug candidate, and 

p then comparing the binding interaction parameter against a predetermined drug 

i'l 

correlation graph (e.g., a mathematical expression fitted to a series of known data 
points) to determine an estimate of at least one pharmacokinetic parameter. Although 

20 many specific details of certain embodiments of the invention are set forth in the 
following detailed description and accompanying Figures, those skilled in the art will 
recognize that the present invention may have additional embodiments, or that the 
invention may be practiced without several of the details described herein. 

For purposes of clarity and to assist in understanding the full scope of the 

25 present invention, a brief review of the nomenclature associated with pharmacokinetics 
has been provided. As used within the context of the present invention, the following 
pharmacokinetic terms shall be construed broadly, and shall have their generally 
accepted meanings as set forth below. (As previously noted, "pharmacokinetics" refers 
to the study of the kinetics associated with the dynamic processes of absorption, 



9 



distribution, metabolism, and excretion (ADME) of a drug and/or its metabolites within 
a living organism.) 

"Absorption" refers to the process of uptake of a drug compound from 
the site of administration into the systemic circulation. The transfer of drug across the 
5 intestinal lumen is generally referred to as oral absorption, whereas the transfer of drug 
across an external physiological barrier is referred to general absorption. As disclosed 
herein, the pharmacokinetic parameter of absorption may be estimated from biosensor 
data associated with a sensor chip having, for example, a plurality of appropriate 
liposomes immobilized thereon. 
10 "Distribution" refers to the transfer of a drug compound from the site of 

administration to the total systemic circulation and then to extracellular and intracellular 
water and tissues. Drug distribution is usually a rapid and reversible process. As 
disclosed herein, the pharmacokinetic parameter of distribution may be estimated from 
biosensor data associated with a sensor chip having, for example, a plurality of 
15 appropriate plasma proteins, liposomes, and/or transport proteins immobilized thereon. 

"Metabolism" refers to the sum of all the chemical reactions for 
biotransformation of endogenous and exogenous substances which take place in the 
living cell. As disclosed herein, the pharmacokinetic parameter of metabolism may be 
estimated from biosensor data associated with a sensor chip having, for example, a 
20 plurality of appropriate metabolic enzymes immobilized thereon. 

"Excretion" refers to the final elimination or loss of a drug from the 
body. Drug excretion includes both passive diffusion and relative specific carrier 
mediated excretion. Drugs may be excreted, unchanged or as metabolites, in urine via 
the kidneys or in feces via the bile and/or the intestine. Volatile compounds are often 
25 excreted in expired air by the lungs. As disclosed herein, the pharmacokinetic 
parameter of excretion may be estimated from biosensor data associated with a sensor 
chip having immobilized thereon, for example, an antibody that specifically detects the 
drug, as well as other proteins/receptors having a high affinity/specificity against the 
drug candidate. Such antibodies and proteins/receptors may be used to quantify the 
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concentration/amount of the drug in different body fluids (e.g., urine/feces) and tissues, 
using a direct binding assay. 

In addition to these ADME parameters, the solubility of a drug is also an 
important property that may be measured by the methods of the present invention. 
"Solubility" refers to the ability of two substances to form a homogeneous solution or 
mixture with each other. Solubility is important for the dissolution of drug given in 
solid dosage form. As disclosed herein, the solubility of a drug candidate may be 
estimated from sensorgram irregularities associated with reflectance minimum and dip- 
shape of biosensor data. 

Furthermore, and as used herein, the term "parameter" refers to a 
constant or variable term in a function (e.g., a mathematical expression) that determines 
the specific form of the function but not necessarily its general nature. For example, the 
constant term "a" in the function f(x) = ax, where "a" determines only the slope of the 
line described by f(x), is referred to as a parameter. As such, the term "binding 
interaction parameter" refers to those constant or variable terms that are related to the 
binding interaction between a drug candidate and a sensing surface-bound biomolecule 
and includes, for example, association and dissociation rate constants, as well as 
maximum binding capacity. Similarly, the term "pharmacokinetic parameter" refers to 
those constant and variable terms that are related to the disposition of the drug candidate 
within a living organism and includes, for example: volume of distribution; total 
clearance; protein binding; tissue binding; metabolic clearance; renal clearance; hepatic 
clearance; biliary clearance; intestinal absorption; bioavailability; relative 
bioavailability; intrinsic clearance; mean residence time; maximum rate of metabolism; 
Michaelis-Menten constant; partitioning coefficients between tissues and blood (or 
plasma) such as those partitioning coefficients associated with the blood brain barrier, 
blood placenta barrier, blood human milk partitioning, blood adipose tissue partitioning, 
and blood muscle partitioning; fraction excreted unchanged in urine; fraction of drug 
systemically converted to metabolites; elimination rate constant; half-life; and secretion 
clearance. 
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The methods of the present invention are intended to be carried out by 
use of an affinity-based biosensor. As is appreciated by those skilled in the art, 
"biosensors" are analytical devices for analyzing minute quantities of sample solution 
having an analyte of interest, wherein the analyte is analyzed by a detection device that 
5 may employ a variety of detection methods. Typically, such methods include, but are 
not limited to, mass detection methods, such as piezoelectric, optical, thermo-optical 
and surface acoustic wave (SAW) device methods, and electrochemical methods, such 
as potentiometric, conductometric, amperometric and capacitance methods. With 
regard to optical detection methods, representative methods include those that detect 

10 mass surface concentration, such as reflection-optical methods, including both internal 
and external reflection methods, angle, wavelength or phase resolved, for example, 
ellipsometry and evanescent wave spectroscopy (EWS), the latter including surface 
plasmon resonance (SPR) spectroscopy, Brewster angle refractometry, critical angle 
refractometry, frustrated total reflection (FTR), evanescent wave ellipsometry, scattered 

15 total internal reflection (STIR), optical wave guide sensors, evanescent wave-based 
imaging such as critical angle resolved imaging, Brewster angle resolved imaging, SPR 
angle resolved imaging, and the like. Further, photometric methods based on, for 
example, evanescent fluorescence (TIRF) and phosphorescence may also be employed, 
as well as waveguide interferometers. 

20 In the detailed description and Examples that follow, the present 

invention is illustrated in the context of SPR spectroscopy. However, it is to be 
understood that the present invention is not limited to this detection method. Rather, 
any affinity-based detection method where an analyte binds to a ligand immobilized on 
a sensing surface may be employed, provided that a change in a property of the sensing 

25 surface is measured and quantitatively indicative of binding of the analyte to the 
immobilized ligand thereon. In the context of SPR spectroscopy, one exemplary class 
of biosensors is sold by Biacore AB (Uppsala, Sweden) under the tradename 
BIACORE® (hereinafter referred to as "the BIACORE instrument"). Such biosensors 
utilize a SPR based mass-sensing technique to provide a "real-time" binding interaction 

30 analysis between a surface bound ligand and an analyte of interest. 
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The BIACORE instrument includes a light emitting diode, a sensor chip 
covered with a thin gold film, an integrated fluid cartridge and a photo detector. 
Incoming light from the diode is reflected in the gold film and detected by the photo 
detector. At a certain angle of incidence ("the SPR angle"), a surface plasmon wave is 
5 set up in the gold layer, which is detected as an intensity loss or "dip" in the reflected 
light. More particularly, and as is appreciated by those skilled in the art, the 
phenomenon of surface plasmon resonance (SPR) associated with the BIACORE 
instrument is dependent on the resonant coupling of light, incident on a thin metal film, 
to oscillations of the conducting electrons, called plasmons, at the metal film surface. 
10 These oscillations give rise to an evanescent field which extends from the surface into 
the sample solution. When resonance occurs, the reflected light intensity drops at a 
sharply defined angle of incidence, the SPR angle, which is dependent on the refractive 
index within the reach of the evanescent field in the proximity of the metal surface. 

Stated somewhat differently, surface plasmon resonance is an optical 
p 15 phenomenon arising in connection with total internal reflection of light at a metal film- 

S liquid interface. Normally, light traveling through an optically denser medium, e.g., a 

^ glass prism, is totally reflected back into the prism when reaching an interface of an 

optically less dense medium, e.g., a buffer, provided that the angle of incidence is larger 
than the critical angle. This is known as total internal reflection. Although the light is 
20 totally reflected, a component of the incident light momentum called the evanescent 
wave penetrates a distance of the order of one wavelength into the buffer. The 
evanescent wave may be used to excite molecules close to the interface. If the light is 
monochromatic and /?-poIarized, and the interface between the media is coated with a 
thin (a fraction of the light wave-length) metal film, the evanescent wave under certain 
25 conditions will interact with free oscillating electrons (plasmons) in the metal film 
surface. When surface plasmon resonance occurs, light energy is lost to the metal film 
and the reflected light intensity is thus decreased. 

The resonance phenomenon will only occur for light incident at a sharply 
defined angle which, when all else is kept constant, is dependent on the refractive index 
30 in the flowing buffer close to the surface. Changes in the refractive index out to about 
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1 from the metal film surface can thus be followed by continuous monitoring of the 
resonance angle. A detection volume is defined by the size of the illuminated area at 
the interface and the penetration depth of the evanescent field. It should be noted that 
no light passes through the detection volume (the optical device on one side of the metal 
5 film detects changes in the refractive index in the medium on the other side). 

As noted above, the SPR angle depends on the refractive index of the 
medium close to the gold layer. In the BIACORE instrument, dextran is typically 
coupled to the gold surface, with the ligand being bound to the surface of the dextran 
layer. (Note a detailed discussion of matrix coatings for biosensor sensing surfaces is 

10 provided in U.S. Patent No. 5,436,161, which is incorporated herein by reference in its 
entirety.) The analyte of interest is injected in solution form onto the sensor surface 
through the fluid cartridge. Because the refractive index in the proximity of the gold 
film depends upon (1) the refractive index of the solution (which is constant) and, (2) 
the amount of material bound to the surface, the binding interaction between the bound 

15 ligand and analyte can be monitored as a function of the change in SPR angle. 

A typical output from the BIACORE instrument is a "sensorgram," 
which is a plot of response (measured in "resonance units" or "RU") as a function of 
time. An increase of 1 ,000 RU corresponds to an increase of mass on the sensor surface 
of approximately 1 ng/mm 2 . As a sample containing the analyte contacts the sensor 

20 surface, the ligand bound to the sensor surface interacts with the analyte in a step 
referred to as "association." This step is indicated on the sensorgram by an increase in 
RU as the sample is initially brought into contact with the sensor surface. Conversely, 
"dissociation" normally occurs when sample flow is replaced by, for example, a buffer 
flow. This step is indicted on the sensorgram by a drop in RU over time as analyte 

25 dissociates from the surface-bound ligand. A detailed discussion of the technical 
aspects of the BIACORE instrument and the phenomenon of SPR may be found in U.S. 
Patent No. 5,3 13,264, which is incorporated herein by reference in its entirety. 

In addition, a detailed discussion of the technical aspects of the biosensor 
sensor chips used in connection with the BIACORE instrument may be found in U.S. 

30 Patent No. 5,492,840, which is incorporated herein by reference in its entirety. This 
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patent discloses, among other things, that each sensor chip may have a plurality of 
sensing surfaces, and that such sensing surfaces may be arranged in series or parallel 
with respect to the fluid sample pathway of the fluid cartridge. This patent also 
discloses that each of the plurality of sensing surfaces of a single sensor chip may have 
5 bound thereto a unique type of ligand that is capable of interacting with an analyte in its 
own characteristic way. 

For example, and as disclosed herein, each of the four discrete sensing 
surfaces of the BIACORE instrument may have immobilized thereon biomolecules such 
as liposomes, plasma proteins, CYP 450 enzymes, other metabolic enzymes, and/or 
10 transport/efflux proteins. By immobilizing one or a selected combination of at least two 
s5 different liposomes, plasma proteins, CYP 450 enzymes, other metabolic enzymes, 

and/or transport/efflux proteins, onto the one or more discrete sensing surfaces of a 
sensor chip, one or more pharmacokinetic parameters associated with a drug candidate 
\U may be readily determined. More specifically, and in one embodiment of the present 

15 invention, it has been discovered that (1) an estimate of an "absorption" parameter of a 
drug candidate may be determined from the binding interactions between the drug 
candidate and one or more appropriate liposomes; (2) an estimate of a "distribution" 
parameter of a drug candidate may be determined from the binding interactions between 
the drug candidate and an appropriate set of plasma proteins {e.g., specific and 
20 nonspecific tissue binding and tissue permeability); (3) an estimate of a "metabolism" 
parameter of a drug candidate may be determined from the binding interactions between 
the drug candidate and an appropriate set of metabolic enzymes; and (4) an estimate of 
an "excretion" parameter of a drug candidate may be determined from the binding 
interactions between the drug candidate and an appropriate set of transport proteins. 
25 Within the context of the present invention, suitable liposomes for 

estimating an "absorption" parameter of a drug candidate include organic compounds 
originating from natural or synthetic lipid molecules such as glycerophospholipids, 
glyceroglycolipids, sphingophospholipids and sphingoglycolipids, and from the classes 
phosphatidyl choline, phosphatidyl etanolamine, phosphatidyl serine, phosphatidyl 
30 glycerol, phosphatidyl acid, phosphatidyl inositol, galactopyranoside, 
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digalactopyranoside, ceramide-phosphatidyl choline, ceramide-phosphatidyl 
etanolamine, ceramide-phosphatidyl serine, ceramide-phosphatidyl glycerol, ceramide- 
phosphatidyl acid, ceramide-phosphatidyl inositol, sphingomyelin molecules, 
glucosylceramides, glucocerebrosides, galactoceramides, galactocerebrosides, 
5 gangliosides, monoacyl phosphatidyl choline, cardiolipin molecules, that may be linked 
to saturated or unsaturated fatty or fluorocarbon chains ranging from eight to twenty- 
four carbons in length where fatty chains attached to the head group can be the same or 
of different structure, cholesterol, lanosterol, ergosterol, stigmasterol, sitosterol and 
derivatives thereof capable of being incorporated into lipid membranes, N,N-dimethyl- 
10 N-octadecyl-l-octadecanammonium chloride or bromide, (N-[l-(2,3 
%3 dioleoyloxy)propyl]-N,N,N-trimethylammonium chloride, diacetyl phosphate, N-[2,3- 

= ij dihexadecy loxy)prop- 1 -y l]-N,N,N^trimethy lammonium chloride, bolaamphiphiles, 

f: Z polyglycerolmonoalkylethers, polyethoxymonoalkylethers, as well as liposome-forming 

*S molecules from the classes amphiphilic polymers, amino acids, crown ether compounds 

s - 15 and di(acyloxy)dialkylsilanes. The liposomes of the present invention and the like may 
h: be immobilized onto the sensing surface by the technique disclosed in Example 1 . 

C3 Suitable plasma proteins for estimating a "distribution" parameter of a 

p drug candidate include proteins such as Immunoglobulin G (7s-y-globulin), IgG; 

™ Immunoglobulin A, IgA; Secretory IgA, s IgA; Immunoglobulin M (19s-y-globulin) 

20 IgM; Immunoglobulin D, IgD; Immunoglobulin E, IgE; a 1 -Antitrypsin, al PI, al A; 
al-Antichymotrypsin, (a 1 X-Glycoprotein) al X; Inter-a-Trypsin inhibitor, laTI; 
Antithrombin III, (heparin cofactor) AT III; ct-Thiol Proteinase inhibitor (LMW 
kininogen) aTPI; C 1 -Inacti vator, CI esterase inhibitor (cc2-neuraminoglycoprotein) CI 
INA: a2-Macroglobulin, a2 M; a2-Antiplasmin, ct2AP; Cystatin C (Post-y-globulin; y- 
25 Trace protein); Clq (1 IS protein); Clr; Cls; C2; C3 (piC-globulin), C4 (plE- 
globulin); C5 (plF-globulin); C6; C7; C8; C9; Factor B, (C3-proacti vator; P2- 
glycoprotein II; glycine-rich P-glycoprotein); Factor D (C3-proactivator convertase); 
Properdin, P; Factor I, (C3b inactivator); C4-binding Protein; Fibrinogen, FI, 
Prothrombin, F II; Factor V (proaccelerin); FV; Factor VII, (proconvertin), F VII; 
30 Factor VIII: C (antihemophilic factor) F VIII: C; Factor VIII - Related Antigen; FV III: 
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Rag (von Willebrand factor) (VWF); Factor IX (Christmas factor) F IX; Factor X, 
(Stuart-Power factor) F X; Factor XI (plasma thromboplalstin antecedent) F XI; Factor 
XII (Hageman factor) F XII; Factor XIII (fibrin stabilizing factor) F XIII; high- 
molecular-weight (HMW) Kininogen (Fitzgerald factor); Prekallikrein (Fletcher factor); 
5 Plasminogen; Protein C; Protein S; Allbumin, ALB; Haptoglobin, HP, Hp 1-1, Hp 2-1, 
Hp 2-2; Prealbumin (transthyretin, thyroxine-binding prealbumin); Retinol-binding 
Protein RBP; Thyroxine-binding Globulin TBG; Transcortin (corticosteroid-binding 
globulin) CBG; Sex hormone-binding globulin (Steroid-binding P Globulin) SHBG; 
Vitamin D-binding Protein (Gc-globulin, group specific component) VDBP; 
r «j 10 Transcobalamin I, TC I; Trancobalamin II, TC II; Transferrin (siderophilin) TF; 
m Ferritin; Hemopexin, HPX; Apolipoprotein A, Apo A-l, Apo A-II; Apolipoprotein B, 

fO Apo B-48, Apo B-100; Apolipoprotein C; Apo C-I, Apo C-II, Apo C-III; 

M- 

Apolipoprotein E, Apo E; Apolipoprotein (a), apo (a); Serum Amyloid A, SAA; ot- 
^ Fetoprotein, AFP; a 1 -Acid Glycoprotein (orosomucoid) al AG; Ceruloplasmin, CP; 

a 15 Serum Amyloid P protein (9.5S a 1 -glycoprotein; al-macroglobulin) SAP; a2-HS 

O 

Glycoprotein, a2 HS; Fibronectin (cold insoluble globulin) FN, C-reactive Protein, 
CRP; p2-microglobulin, p2 M; Pregnancy-specific p 1 -glycoprotein, SP1; al- 
£3 microglobulin, al M The plasma proteins of the present invention and the like may be 

immobilized onto the sensing surface through the use of known immobilization 

20 techniques as is appreciated by those skilled in the art. 

Suitable metabolic enzymes for estimating a "metabolism" parameter 
and/or an "excretion" parameter of a drug candidate include cytochrome P450-enzymes 
selected from the class of cytochrome P450 enzymes having an active site characterized 
by a protoporphyrin IX-iron complex with thiolate from a cysteine of the enzyme 

25 serving as the fifth ligand to iron. Suitable CYP 450 enzymes include CYTOCHROME 
P450, CYP1A1, CYP1A2, CYP2A1, 2A2, 2A3, 2A4, 2A5, 2A6, CYP2B1, 2B2, 2B3, 
2B4, 2B5, 2B6, CYP2C1, 2C2, 2C3, 2C4, 2C5, 2C6, 2C7 ? 2C8, 2C9 ? 2C10 ? 2C11, 
2C12, CYP2D1, 2D2 5 2D3, 2D4 5 2D5 ? 2D6 3 CYP2E1, CYP3A1, 3A2, 3A3 ? 3A4 ? 3A5, 
3A7 ? CYP4A1, 4A2 ? 4A3 5 4A4, CYP4A1 1, CYP P450 (TXAS), CYP P450 11A 

30 (P450scc), CYP P450 17(P45017a), CYP P450 19 (P450arom), CYP P450 51 
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(P45014a), CYP P450 105A1, CYP P450 105B1. Other important metabolic enzymes 
include Glutathione-thioethers, Leukotriene C4,butyrylcholinesterase, human serum 
paraoxonase/arylesterase, N-Acetyltransferase, UDP-glucuronosyltransferase (UDPGT) 
isoenzymes, TL PST, TS PST, drug glucosidation conjugation enzyme, the glutathione- 
5 S-transferases (GSTs) (RX:glutathione-R-transferase), GST1, GST2, GST3, GST4, 
GST5, GST6, alcohol dehydrogenase (ADH), ADH I, ADH II, ADH III, aldehyde 
dehydrogenase (ALDH), cytosolic (ALDH1), mitochondrial (ALDH2), monoamine 
oxidase, MAO: Ec 1.4.3.4, MAOA, MAOB, flavin-containing monoamine oxidase, 
enzyme superoxide dismutase (SOD), Catalase, amidases, Nl,-monogIutathionyl 
fS% 10 spermidine, Nl,N8-bis(glutathionyl) spermidine, Thioesters, GS-SG, GS-S-cysteine, 
^0 GS-S-cysteinylglycine, GS-S-03H, GS-S-CoA, GS-S-proteins, S-carbonic anhydrase 

fU III, S-actin, Mercaptides, GS-Cu(I), GS-Cu(II)-SG, GS-SeH, GS-Se-SG, GS-Zn-R, GS- 

Cr-R, Cholin esterase, lysosomal carboxypeptidase, Calpains, Retinol dehydrogenase, 
z£ Retinyl reductase, acyl-CoA retinol acyltrunderase, folate hydrolases, protein 

s 15 phosphates (pp) 4 st, PP-1, PP-2A, PP-2Bpp-2C, deamidase, carboxyesterase, 

£3 

f2 Endopeptidases, Enterokinase, Neutral endopeptidase E.C.3.4.24.1 1 , Neutral 

t H endopeptidase, carboxypeptidases, dipeptidyl carboxypeptidase, also called peptidyl- 

s s s 

Q dipeptidase A or angiotensin-converting enzyme (ACE) E.C.3.4.15.1, carboxypeptidase 

M, g-Glutamyl transpeptidase E.C.2.3.2.2, Carboxypeptidase P, Folate conjugase 

20 E.C.3.4.12.10, Dipeptidases, Glutathione dipeptidase, Membrane Gly-Leu peptidases, 
Zinc-stable Asp-leu dipeptidase, Enterocytic intracellular peptidases, Amino 
tripeptidase E.C.3.4.1 1.4, Amino dipeptidase E.C. 3.4. 13.2, Prodipeptidase, Arg- 
selective endoproteinase; the family of brush border hydrolases, Endopeptidase-24. 1 1, 
Endopeptidase-2(meprin), Dipeptidyl peptidase IV, Membrane dipeptidase GPI, 

25 Glycosidases, Sucrase-isomaltase, Lactase-glycosyl-ceraminidase, Glucoamylase- 
maltase, Trehalase, Carbohydrase enzymes, alfa-Amylase (pancreatic), Disaccharidases 
(general), Lactase-phhlorizin hydrolase, Mammalian carbohydrases, Glucoamylase, 
Sucrase-isomaltase, Lactase-glycosyl ceramidase, Enzymatic sources of ROM, 
Xanthine oxidase, NADPH oxidase, Amine oxidases, Aldehyde oxidase, 

30 Dihydroorotate dehydrogenase, Peroxidases, Human pancreatic exocrine enzymes, 
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Trypsinogen 1, Trypsinogen 2, Trypsinogen 3, Chymotrypsinogen, proElastase 1, 
proElastase 2, Protease E,Kallikreinogen, procarboxypeptidase Al, 
procarboxypeptidase A2, procarboxypeptidase Bl, procarboxypeptidase B2, 
Glycosidase, Amylase, lipases, Triglycaride lipase, Collipase, Carboxyl ester hydrolase, 
5 Phospholipase A2, Nucleases, Dnase I, Ribonucleotide reductase (RNRs), Wistar rat 
exocrine pancreatic proteins, Label Protein IEP, Al Amylase 1, A2 Amylase 2, Lipase, 
CEL Carboxyl-ester lipase, PL Prophospholipase A, Tl Trypsinogen 1, T2 Trypsinogen 
2, T3 Trypsinogen 3, T4 Trypsinogen 4, CI Chymotrypsinogen 1, C2 
Chymotrypsinogen 2, PE1 Proelastase 1, PE2 Proelastase 2, PC A Procarboxypeptidase 

10 Al, PCA1 Procarboxypeptidase A2, PCB1 Procarboxypeptidase Bl, PCB2 
Procarboxypeptidase B2, R Ribonuclease, LS Lithostatin, Characteristics of UDPGT 
isoenzymes purified from rat liver, 4-nitrophenol UDPGT, 17b-Hydroxysteriod 
UDDPGT, 3-a-Hydroxysteroid UDPGT, Morphine UDPGT, Billirubin UDPGT, 
Billirubin monoglucuronide, Phenol UDPGT, 5-Hydroxytryptamine UDPGT, 

15 Digitoxigenin monodigitoxide UDPGT, 4-Hydroxybiphenyl UDPGT, Oestrone 
UDPGT, Peptidases, Aminopeptidase N, Aminopeptidase A, Aminopeptidase P, 
Dipeptidyl peptidase IV, b-Casomorphin, Angiotensin-converting enzyme, 
Carboxypeptidase P Angiotensin II, Endopeptidase-24. 1 1, Endopeptidase-24. 1 8 
Angiotensin I, Substance P (deamidated), Exopeptidase, 1 . NH2 terminus 

20 Aminopeptidase N (EC 3.4.1 1.2), Aminopeptidase A (EC 3.4.1 1.7), Aminopeptidase P 
(EC 3.4.11.9), Aminopeptidase W (EC 3.4.11.-), Dipeptidyl peptidase IV (EC 3.4.14.5), 
g-Glutamyl transpeptidase (EC 2.3.2.2), 2. COOH terminus Angiotensin-converting 
enzyme (EC 3.4.15.1), Carboxypeptidase P (EC 3.4.17.-), Carboxypeptidase M (EC 
3.4.17. 12),3. Dipeptidase Microsomal dipeptidase (EC 3.4.13.19), Gly-Leu peptidase, 

25 Zinc stable peptidaseJEndopeptidase Endopeptidase-24. 1 1 (EC 3.4.24.11), 
Endopeptidase-2 (EC 3.4.24.18, PABA-peptide hydrolase, Meprin, Endopeptidase-3, 
Endopeptidase (EC 3.4.21.9), GST Al-1, Alpha,GST A2-2 Alpha, GST Mla-la Mu, 
GST Mlb-lb Mu, GST M2-2 Mu, GST M3-3 Mu, GST M4-4 Mu, GST M5-5 Mu, 
GST Pl-1 Pi, GST Tl-1 Theta, GST T2-2 Theta, Microsomal Leukotriene C4 synthase, 

30 UGT isozymes, UGT1.1, UGT1.6, UGT1.7, UGT2.4, UGT2.7, UGT2.11, Pancreatic 
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enzymes, Elastase, Aminopeptidase (dipeptidyl aminopeptidase (IV), Chymotrypsin, 
Trypsin, Carboxypeptidase A, Methyltransferases, O-methyltransferases, N- 
methyltransferases, S-methyltransferases, Catechol-O-methyltransferases, MN- 
methyltransferase, S-sulphotransferases, Mg 2+ -ATPase, Growth factor receptors 
5 Alkaline phosphatase, ATPases, Na, K+ATPase, Ca 2+ -ATPase, Leucine 
aminopeptidase, K + channel. The metabolic enzymes of the present invention and the 
like may be immobilized onto the sensing surface through the use of known 
immobilization techniques as is appreciated by those skilled in the art. 

Suitable transport proteins for estimating an "excretion" (as well as an 

10 efflux, absorption, distribution) parameter of a drug candidate include Glucose a. Na + - 
glucose cotransp GLUT 1, b. Facilitative transporter, GLUT 1-5, GLUT-2, Neutral 
amino acid transporter, Na"-independent system L amino acid transporter, cationic 
amino acid, Y + cationic L-amino acid transporter, Dipeptides FT contransport, 
Nucleosides Na^-dependent and facilitative, Taurine Na" and CI' dependent, Bile acids 

15 Na f /bile acid cotransporter, Na^-independent bile acid transporter, ABC transporters, 
Prostaglandins facilitative transpoprter, NaVFT exchanger Antiporter, Phosphate NaVPi 
cotransporter, Sulfate Na + -cotransporter, Transporters for neurotransmitters, 
Norepinephrine Na7Cl' cotransporter, Dopamine Na7CI~ cotransporter, Serotonine 
Na7Cl* cotransporter, GAB A, GAT-1, NaVCl" dependent, Glycine Na7Cl" dependent, 

20 Glutamate Na + cotransporter, K70H" counter-transport, ABC transporters, P- 
glycoprotein (MDR1, MDR 2 OR MDR 3), CI - channel (CFTR), Antigenic peptides, 
TAP1 and TAP2 heterodimer, Lung resistance protein (LRP), Multidrug resistance 
protein 1 (MRP1), Multidrug resistance protein 2 (MRP2 or cMOAT), Multidrug 
resistance protein 3 (MRP3), Multidrug resistance protein 4 (MRP4), Multidrug 

25 resistance protein 5 (MRP5), Multidrug resistance protein 6 (MRP6), mrp (mouse), 
EBCR (rabbit), C. elegans mrpl (nematode), C. elegans mrp2 (nematode), MRP6 
(human), YCF1 (yeast), AtMRPl (Arabidopsis), SUR1 (human), sur2 (rat, mouse), 
YOR1/YRS1 (yeast), LtpgpA (leishmania), Hepatic amino acid transport system, 
Neutral amino acid transporters, MeAIB, Dicarboxylic amino acids, Neutral amino 

30 acids (branched), b-amino acids, Long chain fatty acids, Monoglycerides, L- 
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Lysophosphatidylcholine, Transport proteins for bilirubin, Bilitranslocase (BTL), 
Organic anion binding protein (OABP), BSP/bilirubin binding protein, Signal receptor 
and transduction-hydrolases, ATPases, Na + dependent/independent bile acid transport, 
Bilirubin/BSP carrier (Cl-dependent), SO/OH exchanger Cl channel Na + /H + exchanger, 
5 NaTHCO 3 cotransport, GSH-, GSSH-, GS conjugate carrier, SOVHCO 3 exchanger, Na + - 
dependent amino acid transport, Dipolar amino acid transporter, Basic amino acids, 
Cystine, Imino acids CI", b- Amino acids Cl\ XAG Acidic amino acids K + , A Dipolar a- 
amino acids, three-and four-carbon dipolar amino acids, L Bulky, hydrophobic, dipolar 
amino acids, y+ Basic amino acids, folate transperter, cbl transport proteins, Na + -K~ - 

10 ATPase, Bile acid transperter (BAT), protein kinase C, NaVP sympoter (NIS); Bile salt 
export pump(BSEP, cBAT, SPGP); Intestinal bile acid transporters; Purine selective 
Na + -dependent nucleoside transporter (hSPNTI); Pyrimidine selective Na + -dependent 
nucleoside transporter (cNTI); Mitoxantrone transporter (MXR1 and MXR2); Intestinal 
oligopeptide transporter(PepTl); Renal oligopeptide transporter(PepT2); Breast cancer 

15 resistance protein (BCRP). The transport/efflux proteins of the present invention and 
the like may be immobilized onto the sensing surface through the use of known 
immobilization techniques as is appreciated by those skilled in the art. 

As stated above, it has been discovered that the pharmacokinetic 
parameters of absorption, distribution, metabolism, and excretion (ADME) of a drug 

20 candidate may be determined from the binding interactions between the drug candidate 
and appropriate sensing surface-bound biomolecules (e.g., from different liposomes, 
plasma proteins, CYP 450 enzymes, other metabolic enzymes, and/or transport/efflux 
proteins, such as those as identified above) of a biosensor. That is, the binding 
interactions between a drug candidate and one or more sensing surface-bound 

25 biomolecules selected from liposomes, plasma proteins, CYP 450 enzymes, other 
metabolic enzymes, and transport/efflux proteins, may be measured via the BIACORE 
instrument to determine a "binding interaction parameter" of the drug candidate. The 
estimated binding interaction parameter may then, in turn, be compared against a 
predetermined drug correlation graph to estimate one or more absorption, distribution, 

30 metabolism, and excretion "pharmacokinetic parameters' 1 of the drug candidate. 
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Moreover, the pharmacokinetic parameter of the drug candidate may be, for example, 
selected from the group consisting of volume of distribution; total clearance; protein 
binding; tissue binding; metabolic clearance; renal clearance; hepatic clearance; biliary 
clearance; intestinal absorption; bioavailability; relative bioavailability; intrinsic 
5 clearance; mean residence time; maximum rate of metabolism; Michaelis-Menten 
constant; partitioning coefficients between tissues and blood (or plasma) such as those 
partitioning coefficients associated with the blood brain barrier, blood placenta barrier, 
blood human milk partitioning, blood adipose tissue partitioning, and blood muscle 
partitioning; fraction excreted unchanged in urine; fraction of drug systemically 
10 converted to metabolites; elimination rate constant; half-life; and secretion clearance; 

*D any one of which may be determined from appropriate sensorgrams of the selected 

.u 

y drug-biomolecule interaction (alternatively, they may be determined by extracting score 

vectors with principal component analysis from digitalized interaction profiles, as well 
as other known multivariate methods). 
15 For example, to estimate an apparent equilibrium constant between a 

drug candidate and selected sensing surface-bound biomolecules, the following 

O procedure may be employed. First, a concentration series {e.g., 0, 20, 50, 100, 500, and 

W 

Q 1,000 |iM) of the drug candidate may be prepared, and sequentially injected into a 

fJ * biosensor having a sensor chip operatively associated therewith, wherein the sensor chip 

20 has a reference sensing surface and at least one sensing surface with surface-bound 
biomolecules. The relative responses at steady-state binding levels for each drug 
concentration level may then be measured. Because of bulk-refractive index 
contributions from solvent additives in the biosensor's running buffer, a correction 
factor may be calculated (via known calibration procedures) and applied to give 
25 corrected relative responses. The corrected relative responses for each drug 
concentration may then be mathematically evaluated as is appreciated by those skilled 
in the art to estimate the apparent equilibrium constant. 

More specifically, the apparent equilibrium constant (KD) between the 
drug candidate and the selected sensing surface-bound biomolecules may be calculated 
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by fitting the measured equilibrium (R cq ) data and known drug concentration (C) data to 
Equation (1): 

Re q = C^JiC + KD) + offset (1) 

5 

wherein is the maximum binding capacity of the sensing surface. Alternatively, 
the apparent equilibrium constant (KD) may be calculated from as little as two 
concentrations (i.e., CI and C2) of the drug candidate by use of Equation (2): 

10 KD = (R eqCI /Cl - R eqC2 /C2)/(R eqC2 - R eqC1 ) (2) 

The estimated apparent equilibrium constant (KD) or binding level at a specific molar 
drug concentration may then, in turn, be compared against a predetermined drug 
fy correlation graph to estimate one or more absorption, distribution, metabolism, and 

p 15 excretion parameters of the drug candidate. 

*0 In the context of the present invention, a predetermined drug correlation 

graph refers to a mathematical expression or function that has been developed from 
binding interaction data associated with known drug compounds. For example, a 
correlation graph may be constructed wherein known binding interaction parameters 
□ 20 (e.g., apparent equilibrium constants and/or binding levels at specific molar drug 
concentrations) for known compounds are plotted along the abscissa (i.e., the x-axis) 
and corresponding measured binding interaction parameters obtained via the biosensor 
are plotted along the ordinate (i.e., the y-axis), or vice versa. Stated somewhat 
differently, the correlation plot is a graphical representation of a mathematical 
25 expression that correlates known and measured affinity data. 

By comparing the estimated binding interaction parameter obtained from 
the selected drug-biomolecule interaction against the mathematical expression (i.e., 
correlation graph) correlated from known and measured affinity data, it has been 
surprisingly discovered that an estimate of a pharmacokinetic ADME parameter related 
30 to the drug candidate may be accurately predicted. In addition, by immobilizing a 
selected combination of at least two different liposomes, plasma proteins, CYP 450 
enzymes, other metabolic enzymes, and/or transport/efflux proteins, onto one or more 
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discrete sensing surfaces of a sensor chip, it has been further surprisingly discovered 
that at least two pharmacokinetic parameters associated with a drug candidate may be 
readily determined. Moreover, by immobilizing a selected combination of different 
liposomes, plasma proteins, CYP 450 enzymes, other metabolic enzymes, and/or 
5 transport/efflux proteins, onto the sensing surfaces of biosensor, such as in a 
predetermined line of different spots of biomolecules, an ADME pattern or profile, such 
as a drug candidate characterization matrix, may be readily developed. Such ADME 
profiles are of great utility for purposes of drug screening. 

In addition to determining one or more pharmacokinetic parameters by 
10 monitoring the refractive index changes of a biosensor as disclosed above, the solubility 

■taF 

%Q of a drug candidate may also be simultaneously determined (together with the one or 

more pharmacokinetic parameters) by monitoring the minimum, maximum or centroid 
of the drug-biomolecule interaction signal (as is disclosed, for example, in PCT 
Publication No. WO 97/09618, which is incorporated herein by reference in its 
15 entirety). More specifically, the solubility of the drug candidate may be 
estimated/identified from irregularities, as well as from reflectance minimum (R™ in ) and 
dip-shape data, shown in the sensorgram of the drug-biomolecule interaction. In 
general, the concentration at which precipitation occurs is referred to as the solubility 
limit; this property may be important to measure because it may indicate that the drug 
20 candidate has an affinity greater than otherwise indicated. 

Solubility problems of the drug candidate may be detected because 
insoluble particles (i.e., precipitates) tend to bind to sensing surface-bound 
biomolecules, thereby causing the sensing surface to be non-homogenous. (A 
homogenous surface has the same surface concentration throughout, whereas a non- 
25 homogenous surface has concentration disruptions.) A non-homogenous sensing 
surface actually measures several refractive indices, which are all averaged together in 
the biosensor's detector. These multiple measurements from a single, but non- 
homogenous sensing surface, tend to result in an increase in reflectance minimum and a 
broadening of the dip associated therewith. 
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An illustration of how sensorgram irregularities may be used to identify 
solubility problems (i.e., the presence of precipitates) is shown in Figures 1A-C. More 
specifically, the steady-state binding levels associated with a selected drug-biomolecule 
interaction is generally shown in Figure 1 A. By enlarging or "zooming" in on the top 
5 area of the sensorgram, which is reflective of the steady-state binding, the sensorgram 
irregularities become more apparent as is shown in corresponding Figure IB. The 
sensorgram irregularities become even more apparent by using reference subtracted data 
(i.e., bulk-refractive index effects have been eliminated) as is shown in corresponding 
Figure 1C. 

™ 10 An illustration of how reflectance minimum (R^J and dip-shape data 

4: may be used to identify solubility problems (i.e., the presence of precipitates) is shown 

fU in Figures 2A-C. More specifically, a homogeneous surface has the same surface 

s > 

lg concentration throughout, and will thus result in a single "dip" with respect to the 

IZ reflected light intensity as shown in Figure 2A. When precipitates bind to the sensing 

if B 

s 15 surface, the sensing surface becomes non-homogenous which tends to result in an 

C3 

m increase in reflectance minimum and a broadening of the dip associated therewith as is 

jr;~ shown in Figures 2B and 2C, respectively. 

C3 Based on the foregoing methods for assaying a drug candidate, 

s 

researchers may now simultaneously measure several different pharmacokinetic 
20 parameters of the drug candidate, as well as gauge the drug candidate's solubility, by 
using a single analytical instrument. The present invention simplifies and improves the 
rate of drug discovery and development because important pharmacokinetic data may 
now be readily obtained at a relatively early stage of the process. 

In addition to the foregoing methods, the present invention is also 
25 directed to apparatuses adapted to carrying out such methods. More specifically, the 
apparatuses of the present invention comprise a biosensor having a sensing surface 
associated therewith, and a computer system that facilitates the implementation of the 
steps associated with the methods disclosed herein. The computer includes a computer 
memory containing a data structure useful for assaying a drug candidate; the data 
30 structure comprises binding interaction data associated with known drug compounds 
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such that the data structure may be used to determine an estimate of at least 
pharmacokinetic parameter of the drug candidate. 

The aspect of the present invention relating to a computer memory 
containing a data structure useful for assaying a drug candidate may be more fully 
5 illustrated in the context of a high-level computer block diagram as is depicted in 
Figure 3. As shown, such a computer system 300 contains one or more central 
processing units (CPUs) 310, input/output devices 320, and the computer memory 
containing a data structure useful for assaying a drug candidate(memory) 330. Among 
the input/output devices is a storage device 321, such as a hard disk drive, and a 
10 computer-readable media drive 322, which may be used to install software products, 
where the software products are provided on a computer-readable medium, such as a 
Hi CD-ROM. The input/output devices also include a network connection 323, through 

9 E which the computer system 300 may communicate with other connected computer 

systems, such as networks. The input/output devices may also contain a display 324 
1 5 and a data input device 325. 

The memory 330 preferably contains an operating system 331, such as 
MICROSOFT WINDOWS, for providing to other programs access to resources of the 
computer system. The memory 330 preferably further contains software 332. While 
the computer memory containing a data structure useful for assaying a drug candidate is 
20 preferably implemented on a computer system configured as described above, those 
skilled in the art will recognize that it may also be implemented on computer systems 
having different configurations. 

In a related aspect, the present invention is also directed to a generated 
data signal conveying a data structure useful for assaying a drug candidate. As above, 
25 the data structure comprises binding interaction data associated with known drug 
compounds, such that the data structure may be used to determine an estimate of at least 
pharmacokinetic parameter of the drug candidate. 

In still a further embodiment, a sensor surface adapted for use with a 
biosensor is disclosed. The sensor surface has a hydrogel matrix coating coupled to the 
30 top surface of the sensor surface, wherein the hydrogel matrix coating has a plurality of 
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functional groups, and at least two different types of liposomes are bonded to the 
plurality of functional groups at discrete and noncontiguous locations on the hydrogel 
matrix coating of the sensor surface (such as disclosed in Example 1). In the context of 
the BIACORE instrument as described above, the sensor surface is preferably in the 
form a sensor chip, wherein the sensor chip has a free electron metal interposed between 
the hydrogel matrix and the top surface of the sensor chip. Suitable free electron metals 
in this regard include copper, silver, aluminum and gold. 

The sensor surface may have a lipophilic substance interposed between 
the different types of liposomes and the plurality of functional groups, wherein the 
lipophilic substance is covalently bonded to the plurality of functional groups. 
Representative lipophilic substances comprise an alkyl chain having from 12 to 24 
carbon atoms, such as stearylamine. Similarly, representative liposomes include 1,2- 
dimyristol-^-glycero-3-phosphocholine (DMPC) and l-palmitoyl-2-oleoyl-s«-glycero- 
3-phosphocholine (POPC). The sensor surfaces of this invention may also have non- 
liposome biomolecules associated therewith. For example, human serum albumin, CYP 
450 enzyme, a metabolic enzyme, and/or a transport protein may be bonded to the 
plurality of may be bonded to the plurality of functional groups at discrete and 
noncontiguous locations on the hydrogel matrix coating of the sensor surface. 

For purposes of illustration and not limitation, the following examples 
more specifically disclose various aspects of the present invention. 

EXAMPLES 
Example 1 

Simultaneous Measurement of Solubility, Plasma Protein Binding, 
llpophilicity and intestinal absorption for three drug candidates a-c 

This example discloses how combined information from each of four 
different flow-cells of a biosensor may be represented in a characterization matrix 
where each of the three drug candidates A-C illustrates a quality pattern {i.e., HSA% 
Bound, Predicted Lipophilicity, Solubility, and Predicted FA%) which is useful for the 
selection of lead drug compounds. 
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Preparation of Sensor Chip 

Three of the four discrete sensing surfaces of a CMS Sensor Chip 
(Biacore AB, Uppsala, Sweden) were modified such that the CMS Sensor Chip had 
surface-bound biomolecules as depicted below in Table 1 . 



TABLE 1 

Surface-Bound Biomolecules of CMS Sensor Chip 



Surface/Cell No. 


Type of Flow-Cell 


Surface Modification 


FC1 


Ref-1 


Unmodified carboxymethyl dextran (CMS) 


FC2 


Target- 1 


Human Serum albumin - HSA (9-12 kRU) 


FC3 


Target-2 


DMPC-hposomes (5-7 kRU) captured on 
stearyl amine 


FC4 


Target-3 


POPC-hposomes (5-7 kRU) captured on 
stearylamine 



More specifically, and to achieve covalent attachment of stearylamine to 
10 surface/cell nos. 3 and 4, the following procedure was employed. A CMS Sensor Chip 
was first inserted into a BIACORE 3000 biosensor (Biacore AB, Uppsala, Sweden). 
The flow of running buffer (isotonic phosphate buffer pH 7.0 (9.6 g Na 0 HPO 4 - 2H 2 0, 
1.7 g KH 2 P0 4 , 4.1 g NaCl to 1 liter)) with 2% dimethylsulfoxide (DMSO) was directed 
to flow-cells 3 and 4 by using appropriate software commands. A mixture of 0.2 M N- 
15 ethyl-N , -(3-diethylaminopropyl)-carbodiimide and 50 mM N-hydroxysuccinimide in 
water was then injected into the biosensor so as to flow over flow-cells 3 and 4 for a 
period of 10 minutes. The CMS Sensor Chip was then washed with running buffer, 
removed from the biosensor, and placed in a petri dish on top of a few layers of tissue 
paper that had been moistened with ethanol. The exposed sensing surfaces were then 
20 treated with 30 ^il of a 10 mM stearylamine solution in 99% ethanol. After another 45 
minute period, the CMS sensor chip was gently washed with ethanol and then with 
water. Next, the exposed sensing surfaces were further treated with 50 jal of a 1M 
ethanolamine solution at pH 8.5. Finally, and after another 10 minute period, the 
sensing surfaces were washed with HBS-buffer. 
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Following covalent attachment of stearylamine to surface/cell nos. 3 and 
4 - and after the CMS Sensor Chip that had just been washed with HBS-buffer was 
blown dry with nitrogen and reinserted into the biosensor - human serum albumin 
(HSA), DMPC-liposomes, and POPC-liposomes (available from Sigma) were then 
captured onto surface/cell nos. 2, 3 and 4, respectively. The following procedures were 
employed to capture these biomolecules. 

To capture human serum albumin the flow was initially directed to flow- 
cell 2 by using appropriate software commands. The sensing surface of flow-cell 2 was 
then activated for a period of 7 minutes with EDC/NHS. Next, human serum albumin 
was injected at 15 |ig/ml in lOmM acetate buffer pH 5.2 for a period of 7 minutes; then 
1M ethanolamine pH 8.5 was injected for another period of 7 minutes. This procedure 
resulted in the immobilization of 9-12 kRU of human serum albumin to the sensing 
surface of surface/cell no. 2. 

To capture the two different liposomes the flow was first directed to 
flow-cell 3, and 0.5-1 mM DMPC liposome was injected until 5-7 kRU of the DMPC 
liposome had been captured. The flow was then directed to flow-cell 4, and 0.5-1 mM 
POPC liposome was injected until 5-7 kRU of the POPC liposome had been captured. 
Finally, the flow-system of the biosensor was washed with lOOmM NaOH, and the 
autoinjector tubing was washed with 0.5% SDS and 50 nM glycine pH 9.5 to remove 
trace lipids. 

Calibration Procedure 

In order to improve data quality by reducing bulk-refractive index 
contributions from the DMSO solvent additive, a calibration procedure was employed. 
(Note the calibration procedure required varying concentrations of DMSO.) Because 
the running buffer contained 2% DMSO, a series of 8 to 10 calibration solutions were 
made that had varying DMSO concentrations ranging from 1.5% to 3.0%. Each 
calibration solution was sequentially injected over each of the four sensing surfaces by 
using the serial injection mode of the biosensor, and the respective steady-state binding 
levels were measured. The relative responses of the reference sensing surface of flow- 
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cell 1 (FC1 having unmodified CMS as a reference surface) were then subtracted from 
the corresponding calibration responses of the target sensing surfaces of flow-cells 2, 3, 
and 4, respectively, and plotted as functions of the relative responses. From these plots, 
calibration functions were calculated for each of the respective sensing surfaces. 
5 The calibration functions were then used to calculate appropriate 

correction factors for the samples containing the three drug candidates A-C. That is ? for 
the samples containing the drug candidates, the relative responses of the steady state 
binding levels of the reference sensing surface (i.e., signals from flow-cell 1) were 
measured and the respective calibration functions were used to calculate correction 
10 factors appropriate for each drug candidate. The correction factors were then applied 
(i.e., subtracted) to the differences between the reference responses and the sample 
responses to give a responses where the bulk-refractive index contributions from the 
DMSO solvent additive had been eliminated. 

Corrected Relative Responses of Drug Candidates A-C 

15 The modified CMS Sensor Chip (having surface-bound biomolecules as 

depicted above in Table 1) was used to measure the binding interaction of the three drug 
candidates A-C. Concentration series (i.e., 0, 20, 50, 100, 500, and 1,000, \xM) of each 
drug candidate were prepared, and sequentially injected in serial mode into the 
biosensor such that four sensorgrams corresponding to each flow-cell were developed. 

20 Corrected relative responses at steady-state binding levels for each drug concentration 
were determined for each of the target sensing surfaces/flow-cells; such data is 
presented below in Tables 2-4. 
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TABLE 2 

Drug Concentration Data for Human Serum Albumin 



Candidate 


10 nM 


50 nM 


100 ^iM 


500 nM 


1,000 nM 


A 


50 


110 


130 


140 


145 


B 


0 


0 


0 


10 


50 


C 


10 


50 


100 


110 


115 
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TABLE 3 

Drug Concentration Data for POPC Liposome 



Candidate 


10 nM 


50 [iM 


100 |iM 


500 \iM 


1,000 nM 


A 


10 


300 


400 


450 


500 


B 


0 


0 


40 


60 


70 


C 


0 


0 


0 


0 


10 



5 TABLE 4 

Drug Concentration Data for DMPC Liposome 



Candidate 


10 


50 ]iM 


100 ^iM 


500 |iM 


1,000 nM 


A 


0 


50 


70 


90 


95 


B 


0 


10 


20 


50 


55 


C 


0 


0 


0 


0 


0 



10 Analysis of Data 

The corrected relative responses for each drug concentration as depicted 
above in Tables 2-4 were, among other things, adjusted for differences in molecular 
weight of the drug candidates, as well as transformed to increase correlation with other 
properties. For example, for each target sensing surface/flow-cell, an apparent 

15 equilibrium constant (KD) was calculated by fitting the measured equilibrium (R^) data 
and known drug concentration (C) data to Equation (3): 

= C*R max /(C + KD) + offset (3) 

20 wherein is the maximum binding capacity of the sensing surface. In the fitting 
procedure, R^, KD, and the offset were calculated. 

Alternatively, the apparent equilibrium constant (KD) may have been 
calculated from two concentrations of the drug candidate, CI and C2, by use of 
Equation (4): 

25 
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KD = (R eqC1 /C 1 - Re qC2 /C2)/(R eqC2 - R eqCI ) ( 4 ) 

In either case, the KD values obtained via the biosensor were 
subsequently correlated with known KD values calculated from reported human serum 
albumin binding percentages (i.e., reported binding percentages that have been reported 
for known drug compounds), thereby enabling a prediction of the degree of protein 
binding for each of the three drug candidates A-C. In other words, a correlation graph 
as shown in Figure 4 was constructed having known KD values for known compounds 
plotted along the abscissa (i.e., the x-axis) and corresponding measured KD values 
obtained via the biosensor plotted along the ordinate (i.e., the y-axis). This correlation 
graph was then used to predict the degree of plasma protein binding for each of the 
three drug candidates A-C (see below). 

In addition, for each target sensing surface/flow-cell, a molecular weight 
adjusted response at a single concentration level was used as a threshold value for 
ranking the drug candidates A-C. More specifically, a correlation graph as shown in 
Figure 5 was constructed, wherein respective binding levels at 100 (R 100 uM) 
divided by the known drug compounds' molecular weight were plotted along the 
abscissa (i.e., the x-axis) and corresponding human serum albumin binding percentage, 
as measured by equilibrium dialysis, were plotted along the ordinate (i.e., the y-axis). 
This correlation graph was then used to discriminate between strong and weak human 
serum albumin binders for each of the three drug candidates A-C. 

Furthermore, for each target sensing surface/flow-cell, reference 
subtracted sensorgram traces (e.g., FC2 minus FC1) were developed to detect the 
presence of precipitates (caused by low solubility). More specifically, reference 
subtracted sensorgram traces for each of the three drug candidates A-C, as shown in 
Figure 6, were constructed. (Note that when precipitates form and then bind/associate 
with the sensing surface, the response signal will decrease in a non-continuous way.) 
By analyzing the concentration series for each of the three drug candidates A-C (as 
shown in Figure 6), it was determined that drug candidate A precipitated and thus had 
relatively low solubility, whereas drug candidates B-C did not precipitate and thus had 
relatively high solubilities. 
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Finally, and as shown below in Table 5, the combined information from 
each target sensing surface/flow cell was tabulated into a reduced characterization 
matrix where each drug candidate A-C received a quality pattern that was useful for 
selection of the lead drug compounds. (Note that the predicted lipophilicity and fraction 
absorbed are prophetic.) 



TABLE 5 

Characterization Matrix for Drug Candidates A-C 



Candidate 


HSA % Bound 


Pred. Lipophilicity 


Solubility 


Predicted FA% 


A 


99.9 


4 


<1 uM 


>97 


B 


<90 


2 


OK 


>90 


C 


>90 


-2 


OK 


<50 



Based on the forgoing characterization matrix, it was determined that 
drug candidate B was the preferred compound due to its low plasma protein binding 
level (HSA %bound <90), medium lipophilicity (predicted lipophilicity = 2), no 
identified solubility problems (solubility = OK), and acceptable intestinal absorption 
(predicted fraction absorbed >90). 



EXAMPLE 2 

Demonstrated Correlation Between Biosensor Data and 
Fraction Absorbed in Humans (FA%) 

This example discloses a correlation between biosensor data obtained 
form the BIACORE instrument (i.e., BIACORE 3000) and known data for the fraction 
absorbed in humans (FA%) for a number of different drugs, wherein the correlation 
graph is useful for drug candidate absorption predictions. More specifically, a 
correlation graph as shown in Figure 7 was constructed having known fraction absorbed 
in humans (FA%) plotted along the ordinate (i.e., the y-axis) and corresponding 
calibrated (i.e., reference subtracted) steady state binding levels for each drug at 500 
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|iM plotted along the abscissa in a 10-logarithm scale {i.e., the x-axis). In this example, 
the sensing surfaces of the target flow-cells each had 6,000 RU of POPC-GM3 
ganglioside (available from Sigma) captured on stearylamine tiles. (Note that an 
unmodified sensing surface of a CMS Sensor Chip was used as the reference.) As 
5 shown in Figure 7, there is high degree of correlation between the biosensor data and 
the known fraction absorbed in humans (FA%) data, as is evidenced by the 
mathematical expression/function that has been fitted to the various data points. 

Moreover, the correlation graph also shows a classification of the various 
substances, wherein absorption of the substances using the passively transported trans- 

10 cellular route through the intestine are depicted by h = high, m = medium, and 1 = low, 
and wherein absorption using active transport is depicted by h-act = high, and wherein 
absorption of substances with molecular weights <200 using the para-cellular route is 
depicted by h-para = high and m-para = medium, respectively. (Note the Sulfasalazine 
is a pre-drug which rapidly decomposes in the intestine, which may explain its outlier 

15 properties; Piroxicam has a very low solubility ««500 ^M, which may explain its 
lack of correlation.) 

While the present invention has been described in the context of the 
embodiments illustrated and described herein, the invention may be embodied in other 

20 specific ways or in other specific forms without departing from its spirit or essential 
characteristics. Therefore, the described embodiments are to be considered in all 
respects as illustrative and not restrictive. The scope of the invention is, therefore, 
indicated by the appended claims rather than by the foregoing description, and all 
changes that come within the meaning and range of equivalency of the claims are to be 

25 embraced within their scope. 
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